Solitary and compact-like shear waves in the bulk of solids 
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We show that a model proposed by Rubin, Rosenau, and Gottlieb [J. Appl. Phys. 77 (1995) 
4054], for dispersion caused by an inherent material characteristic length, belongs to the class of 
simple materials. Therefore, it is possible to generalize the idea of Rubin, Rosenau, and Gottlieb to 
include a wide range of material models, from nonlinear elasticity to turbulence. Using this insight, 
we are able to fine-tune nonlinear and dispersive effects in the theory of nonlinear elasticity in order 
to generate pulse solitary waves and also bulk travelling waves with compact support. 
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The emergence of significant solitary waves from the 
equations governing the motion of solids requires usually 
a fine balance between nonlinearity and dispersion. The 
first ingredient finds its natural source in the theory of 
finite, fully nonlinear, elasticity, but the source for the 
second is more elusive. Often, specific boundary condi- 
tions (such as: a thin film coated on a substrate, the finite 
interface between phases in an elastic composite material, 
the lateral dimensions of a rod with free surfaces, etc.) 
are used because they introduce a characteristic length, 
to be compared with the wavelength. For homogeneous 
infinite (bulk) solids, the presence of a microstructure is 
invoked, recording the presence of inherent material char- 
acteristic lengths due to, for instance: the atomic lattice 
cells in a metal, the grain size in an alloy, the persistence 
and contour length in a polymer, etc. (it is worth not- 
ing that a similar situation also occurs in fluid dynamics, 
when inherent characteristic lengths may be not ignored 
in models of turbulence.) Although the standard consti- 
tutive equations of continuum mechanics can accommo- 
date large strains and nonlinear material responses, they 
do not in general include the effects of characteristic ma- 
terial lengths 0. The introduction of such information 
in models of material behavior is a subtle and difficult 
problem 0], and the modifications of the basic constitu- 
tive equations proposed in the literature are often unsat- 
isfactory because they lead to physically unrealistic and 
mathematically ill-posed models Q, Q • Finally, finding 
a simple connection between nonlinear and dispersive ef- 
fects in material modeling is an interesting problem per 
se, because the balance between these two effects is at the 
origin of many fascinating phenomena, and in particular 
the propagation of solitons and compactons |^ . 

In 1987 Rosenau [5j proposed a regularized functional 
expansion of the equations for the dynamics of dense 
lattices in order to yield the leading effect of disper- 
sion due to a characteristic length, via a physically rea- 



sonable and mathematically tractable differential equa- 
tion. One of the advantages of this equation over other 
micro-mechanical or nonlocal approaches is that it mod- 
els dispersion without increasing the number of initial 
or boundary conditions of the usual wave equation. On 
the same pathway, a three-dimensional phenomenological 
continuum model of the Rosenau approach was proposed 
by Rubin, Rosenau, and Gottlieb (RRG) in 1995 3]. One 
aim of this Letter is to put in perspective that latter 
model by showing that it is a special case of the theory 
of simple materials in the sense of Noll 0, hence gen- 
eralizing the original idea of Rosenau to a large class of 
material behavior, from nonlinear elasticity to fluid tur- 
bulence. Another objective is to find a clear and direct 
way to balance nonlinearity and dispersion in order to 
generate solitary waves with compact support in solids. 
This balance is attained by relying solely on constitutive 
arguments, without any resort to reductive asymptotic 
expansions; it is then compared (favorably, as it turns 
out) with the results obtained by such methods. 

Let the motion of a body be described by a relation 
x = x(X,t), where x denotes the current coordinates of 
a point occupied by the particle of coordinates X in the 
reference configuration at time t. The deformation gra- 
dient F(X,t) and the spatial velocity gradient L(X,t) 
associated with the motion are defined by F = Grad X 
and L = grad v (v = dx/dt), respectively (in this Let- 
ter, grad and div denote respectively the gradient and 
divergence operators with respect to the current position 
x, and Grad and Div denote the corresponding opera- 
tors with respect to X.) The mathematical theory of 
simple materials is well-grounded from the mathemati- 
cal and thermomechanical points of view |(|; there, the 
Cauchy stress T is determined by the whole history of 
the deformation gradient, 



T = GZo (F(t - a)) . 



(1) 
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where Q is the constitutive functional. We now summa- 
rize the RRG model, which itself follows the axiomatic 
approach to continuum mechanics introduced by Green 
and Naghdi 0. Let us introduce the specific entropy 
per unit mass 77, the specific rate of internal produc- 
tion of entropy £, the entropy flux per unit of present 
area p, and the specific Helmholtz free energy ip. The 
set {r), ip, T} must satisfy, for all possible thermome- 
chanical process, the restrictions coming from the bal- 
ance of angular momentum, the energy equation, and 
from a suitable statement of the second law of thermo- 
dynamics. The central idea of Rubin et al. 3] is to 
add terms {ip 2 ,T 2 }, modeling material dispersion, to the 
Cauchy stress tensor and the free energy: T = T\ + T 2 , 
ip = tpi +^2: say, while the other quantities p} and 
the associated constitutive equations stay unchanged. 
They show that the balance of angular momentum and 
the energy equation then lead to T 2 



T 2 T and 



piP 2 -T 2 -D = 0, 



(2) 



respectively, where p is the mass density per unit volume, 
a superposed dot denotes the material time derivative (at 
X fixed), and D = (L + L T )/2 is the stretching tensor. 
Further, the RRG model considers that ip 2 is a single- 
valued function: tp 2 — ip2(S), where S = D • D is an 
isotropic invariant; then from follows that 

T 2 = pip' 2 (5) [grad v + (grad v) T + 2L T L - 4£) 2 ] . (3) 

Our first result is to show that T 2 is a special case of the 
Cauchy stress er associated with those non-Newtonian 
simple fluids denoted in the literature as second-grade 
fluids 



cr = vAi + a\A 2 + a 2 A\, 



(4) 



where v is the classical Newtonian viscosity, a\, a 2 are 
the microstructural coefficients (which may be taken as 
constants or as functions, for example of 5), and the first 
two Rivlin-Ericksen tensors A%, A 2 are defined by: A\ = 
2D, A 2 = Ai + AiL + L T A 1 . Indeed, recall that L — 
(Grad v)F~ x = FF^ 1 , from which follows that 

L = (Grad u)^ 1 + £F(F Ll ) = grad v - L 2 , (5) 

(where the last equality comes from the material time 
derivative of the identity FF^ 1 = 1). Substituting into 
Ai = L + L T and into the definition of A 2 , we find 
that the tensors J3J and (J2J coincide when v = and 
oti(> 0) = — a 2 = pip' 2 - We conclude that when T\ is in 
the material class then the total Cauchy stress tensor 
T = T\ + T 2 of the corresponding RRG model is also 
that of a simple material. This observation has many in- 
teresting consequences. For example, if p is the Lagrange 
multiplier associated with the constraint of incompress- 
ibility (det F = l,tr D = 0), then the classical choice 
of the Navier-Stokes stress tensor T\ = —pi + vD for 



the RRG model leads to those incompressible homoge- 
neous fluids of second grade for which a± + a 2 = 0. We 
note that the equations of motion associated with those 
fluids, in the general and some special case, have been 
found to coincide with the equations derived by other 
means (asymptotic expansions, Lagrangian means, aver- 
age Euler equation, etc.) in several models of turbulence 
and shallow water theory [JJj- It is now clear that the 
presence of the RRG characteristic length is the reason 
of such connections. 

Let us consider the case of the RRG model when T\ 
corresponds to a hyperplastic, incompressible, isotropic 
solid 



Ti = -pi + 2(dZ/dh)B - 2(dZ/dI 2 )B- 



(6) 



Here B = FF T ', and S = H{I\,I 2 ) is the strain en- 
ergy density, a function of the invariants I\ — tr B and 
I 2 = tr B 1 . Then the RRG model is a special case of 
the so-called solid of second grade [HI . A rectilinear 
shearing motion is defined by 



x = [X + u(Y, t)} ei + Ye 2 + Ze 3 



(7) 



where (ei, e 2 , 63) is a fixed orthonormal triad. Then B — 
l + tiyei®ei + uy(ei®e 2 + e 2 (g)ei), B 1 = 1 + Uye 2 ® 
e 2 - u Y (ei ® e 2 + e 2 ® ei), A\ = uy t (ei (&e 2 + e 2 ® e%), 
and A 2 = u Y tt(ei ® e 2 + e 2 ® ei) + 2u\ t e 2 ® e 2 . Also, 
I\ = I 2 = 3 + u Y and S — u Yt - The equations of motion 
in the absence of body forces (div T — pv) reduce to the 
three scalar equations [l2|: 



dp 
dX 



dY 



d 2 u dT 22 
' dt 2 ' dY 





dZ 



0, (8) 



for which we used ©, or its equivalent form: T 2 = 
pi// 2 (A 2 — Af). The third equation gives p — p(X,Y,t) 
and then the cross-differentiation of the first and sec- 
ond equations leads to p(X, Y, t) — —2Y> 2 u Y + pip 2 u Yt + 
Pi(t)X + po(t), where po, p\ are arbitrary functions of t 
alone and X, is the partial derivative of X with respect to 
Ii. We are left with a single determining equation for the 
displacement u(Y, t): p\ + [2(Ei + E 2 )uf + pip' 2 UYtt)Y — 
putt or, taking the derivative with respect to Y , a single 
equation for the strain U = uy, 



pU + 0U tt Yy = pU tu A = 2(Ei + E 2 ), = prp' 2 . (9) 



Here p — p{U 2 ) > is the generalized shear modulus of 
nonlinear elasticity [l3lll^ |: it is a constant (Lame's sec- 
ond constant) in linear elasticity, and also in nonlinear 
elasticity for the special cases of the neo-Hookean ("ge- 
ometric" nonlinearity) and of the Mooney-Rivlin ( "geo- 
metric" and "physical" nonlincarities) forms of the strain 
energy W (= Ci (/i - 3), Ci (Ji - 3) + C 01 (I 2 - 3), re- 
spectively). However in general, a nonlinear constitutive 
equation for the solid yields a non-constant p; similarly, 



FIG. 1: Profile of the solitary wave when pc 2 /po = 1.1 (thick 
curve), 2 (thin curve), 3 (dashed curve). 



FIG. 2: Profile of the compact wave. 



a general dispersive material behavior is attested by a 
non-constant function = (3{Uf) > (the case where p 
and (3 are both constants is considered elsewhere [l2|). 
To illustrate some features of the RRG model, we make 
the following constitutive choices: p{U 2 ) = po + HiU 2 
and a(U 2 ) — p(/3o + PiUf), where > is the infinites- 
imal shear modulus of the solid. We seek a travelling 
wave solution: U = V(£), £ = Y — ct, with drop bound- 
ary conditions: V r( "- ) (±oo) = 0, n > 0. Then successive 
integrations of the reduced governing equation lead to 

pc 4 /3i(F') 4 + 2pc 2 (3 (V') 2 + [lhV 2 - 2{pc 2 - p )]V 2 = 0. 

(10) 

For a linear variation of the dispersion free energy tp2 
with the isotropic invariant 5 = D ■ D, we have /3i = 
and the following solitary wave solution: 



V(H) 



per - p 



sech 



Mi 



I PC 2 - Lip £ - £ 

pc 2 VFo ' 



(11) 



where £o is arbitrary. This wave is supersonic with re- 
spect to the speed of an infinitesimal shear wave but c is 
otherwise arbitrary. Figure 1 shows its profile for several 
values of pc 2 /po. 

Now, we know from previous results 0] that (|10|) has 
no compact-like solutions when flofli ^ 0. On the other 
hand, when /3q = 0, /3i ^ 0, it is of the form 



(V r f = (70 - 7i V 2 )V 2 



(12) 



(where 70 = 2(pc 2 - p ) / (pc 4 /3i) and 71 = pi / '(pc 4 , 
the right hand-side of which has a double zero (at V = 0) 
and two simple zeros (at V = ± ^70 /71 ) ■ This ob- 
servation alone is sufficient to conclude that (|12fl ad- 
mits a compact wave solution. As it happens, l|12|l can 
be integrated in terms of special functions. Indeed, in- 
troducing the hypergeometric function Hyp 2 Fi[a, b, c, z], 
solution to the second-order linear differential equation 
z(l - z)y" + [c - (a + b + l)z]y' - aby = 0, we find the 
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following identity, 

J (u 2 - u^y^Wu = 2(ii 2 ) 1 / 4 Hyp 2 F 1 

'(13) 

We call l{u) that function and observe that it is defined 
on the interval < u < n/ \[2, where it increases in a 
monotone way from to 1. Then we manipulate l|12|) to 
give 



V / 7i77oV(£) 



(14) 



where £0 is arbitrary. Hence we can build a weak solution 
with finite support measure \/2ttj 1 , defined by: 

Vt77toV(£) = ill' 4 on [0, tt/o], 

= I- 1 [7i 1/4 (V2vr-0] on [7r/a,V27r/a], 

(15) 

and zero everywhere else, see Figure 2. (Here a — 
(471 ) 1 / 4 ). This wave shares some characteristics with 
the solitary wave (|llf) : it is supersonic, its amplitude 
is proportional to y/^o/li = V 2 (P c2 _ Mo)/Mi, growing 
with the speed; in contrast, its width evolves as r y 1 1 ^ 4 , 
which is proportional to c. We emphasize that this com- 
pact wave is derived solely on the basis of constitutive 
arguments, by taking the generalized shear modulus p 
linear in the squared shear (nonlinear effect) and the free 
energy density tp2 proportional to 5 2 (dispersive effect). 
Previous attempts at finding compact waves relied on ad 
hoc asymptotic reductive approaches 0] and to the best 
of our knowledge, a bulk compact wave has never been 
generated in the framework of nonlinear elasticity 0, p. 
91] 

Finally we look at reductive asymptotic expansions 
of the strain equation of motion ©. We consider that 
the generalized shear modulus is slightly nonlinear, p = 
Po(l + ell 2 ) where e is a small parameter, and that the 
dispersion is also of order e: for instance, when ip2 is 
linear in S, we take: (3 = 2(pa 2 )e, where a is a con- 
stant with the dimensions of a length. With the stan- 



4 



dard semi-characteristic variable stretching transforma- 
tion [13: U(Y,t) = vfaer), where r = [^ /(/oa 2 )] 1/2 i 
and rj = Y/a — r are dimensionless variables, we find 
that at order e, © is the modified KdV equation, 

V T + ^(V 3 ) V + V nm = 0, (16) 

for which the solitonic solution is a sech- wave as in (|11H . 
When ip2 is proportional to 8 2 , we take (3 in the form 
/3 = 2(p 2 a 4 /p Q )eU 2 , and find 

«r + 3(A + sK)<W=°. 

an equation similar to the K(3, 3) KdV equation [3.Hj|. 
Although Ijl6(l and (|17|l might have rich implications, we 
do not study them further because the generality of 
has clearly been lost in the process of deriving them. In 
particular, the traveling equation (|l(Jfl is just as reward- 
ing, and is obtained exactly, without the need for a time- 
stretching variable, a moving frame, and an asymptotic 
expansion where higher-order terms are simply ignored. 

In conclusion, we believe that this investigation demon- 
strates the usefulness, the versatility, and ultimately the 
beauty of the RRG model. Although the RRG model 
originates from a microstructural model |5j , it can now be 
apprehended as a purely phenomenological approach to 
dispersion, applicable to the whole class of simple materi- 
als. It overcomes the usual problems of nonlocal theories, 
it unifies several results in the literature, and it provides 
a natural and detailed study of localization of traveling 
waves in elastic materials. Yet several topics raised here 
merit further scrutiny: is the solitary wave exhibited a 
soliton? is the compact wave a compacton? Also, the 
analysis relied on specific constitutive equations. Hence 
for the nonlinear elasticity, the generalized shear modu- 
lus was that of a solid having strain energy density in the 
form E = Cio Ji + C* i J 2 + C 20 Jf + C u Ji J 2 + C Q2 Jf 
where J\ = I\ — 3 and J 2 = I 2 — 3 (so that po = 
2(Ci + Cox) and p x = 4(C 20 + C n + C 02 )). Although 
this, and higher-order, truncation of the polynomial ex- 
pansion of £ are highly popular in experimental nonlinear 
elasticity, in curve fitting, and in finite element packages, 
it should be acknowledged that the determination of the 
"best" Cij is a messy procedure |20J . Eventually, it might 
be of greater (and safer) interest to consider other types 
of nonlinear strain energy densities, such as those aris- 
ing in biological soft tissues, which incorporate the effect 
of strain hardening. Similar remarks hold for the non- 



linear dispersion, for which the free energy was taken as 
ijj 2 = a 5 + (a 1 /2)5 2 . 
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